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Abstract 

A class of dynamical models of turbulence living on a one-dimensional 
dyadic-tree structure is introduced and studied. The models are ob- 
tained as a natural generalization of the popular GOY shell model of 
turbulence. These models are found to be chaotic and intermittent. 
They represent the first example of (l+l)-dimensional dynamical sys- 
tems possessing non trivial multifractal properties. The dyadic structure 
allows to study spatial and temporal fluctuations. Energy dissipation 
statistics and its scaling properties are studied. Refined Kolmogorov 
Hypothesis is found to hold. 



2 On leave on absence from Dipartimento di Fisica, Universita di Tor Vergata. 
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1 Introduction 



Spatio-temporal intermittency is the most intriguing aspect of a fully devel- 
oped three-dimensional turbulent flow. Turbulent structures are thought to be 
generated by chaotic intermittent energy transfer from large to small scales. 
The cascade is pictorially described by the Richardson scenario: large scale 
eddies destabilize and generate small scale eddies with shorter eddy turnover 
times. In this way, a hierarchy of fluctuations on smaller and smaller scales 
and with shorter and shorter characteristic times is produced. 
Kolmogorov 1941 theory |lj describes the statistics of velocity differences, 5 r v, 
at scale r in terms of the averaged energy dissipation e, neglecting completely 
possible spatio-temporal fluctuations. The velocity field statistics is charac- 
terized, among others, by the scaling exponents, C(p)> °f structure functions 
S p (r), in the inertial range: 

S p (r) =< \v(x + r) - v{x)\ p >=< \5 r v\ p >~ r c(p) . (1) 

In the Kolmogorov description, a simple dimensional argument leads to the 
predictions: 

< \5 r v\ p >~ £ p/V/ 3 ; C(P) = P/3. 

On the other hand, experiments show |2j that scaling exponents deviate from 
the linear behaviour. This departure from linearity is the main signature of 
intermittency and implies non-gaussian probability distribution functions for 
the velocity differences in the inertial range. 
Intermittency also affects energy dissipation statistics. 

Experiments || show that the energy dissipation defines a multifractal measure 
on the fluid volume. The multifractal measure is characterized by the scaling 
properties of the coarse-grained energy dissipation, e r , namely: 

d 3 xe(x)] \ ~ r r(p \ (2) 

A(r) J I 

where < ■ > means averaging over all boxes A(r) of size r in which the volume 
occupied by the fluid is partitioned. 

The celebrated Refined Kolmogorov Hypothesis (RKH) [Q links the statistics 
of velocity differences in the inertial range with the statistics of the coarse- 
grained energy dissipation: 

~ e r (6) 

r 

where the symbol ~ means that quantities on both sides have the same scaling 
properties. 
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Using (^|), it is easy to relate the scaling exponents £(p) with the scaling expo- 
nents of the energy dissipation: 

C(p)=r(p/3)+p/3. (4) 

This relation is the natural consequence of the RKH and it is quite well satisfied 
experimentally ||. On the other hand, no satisfactory theoretical arguments 
which predict (H) have ever been found. 

The simplest way to explain phenomenologically the presence of intermittent 
deviations consists in describing the energy transfer mechanism in terms of 
fragmentation stochastic processes. In these models || |6], [7|, one introduces a 
set of eddies leaving on a dyadic structure. 

Random fragmentation models state that the energies contained in eddies at 
successive scales are connected eachother by independent stochastic variables 
with distribution that does not depend on the scale. For example, in the 
random-/3-model we have e n = /3s n+ i, where e n denotes the energy transfer due 
to a typical eddy of size r n = 2~ n r . In the random-/3-model also the active 
volume, V n , occupied by eddies of size r n , is supposed to change randomly from 
scale to scale: V n+ i = j3V n . From the definition (g) we have r{p) = — log 2 < 
j3' p ~ l > and the corresponding expression £(p) = p/3 — log 2 < j3 p l z ~ l > must be 
valid if the RKH holds. Due to the freedom in the choice of the f3 distribution 
all these models are able to fit with good accuracy experimental values. On 
the other hand, they lack any direct linking with the original Navier-Stokes 
(NS) equations. 

In this paper, we investigate a class of dynamical models which fill the gap 
between purely stochastic fragmentation models and the original NS dynamics. 
In particular, we define and study a dynamical model on a (l + l)-dimensional 
(that is, having space and time dimensions) dyadic tree. We decompose the 
original velocity field in terms of fluctuations localized both in Fourier and 
real spaces. One can look at this model as an approximation of the original 
NS equations in a wavelets basis (see next section). Several similar models 
have been already investigated in the past by Grossmann and coworkers |8]. [5j . 
At difference from models studied in || |, we do not impose all the original 
geometrical constraints of the NS equations. For example, we do not have 
neither divergenceless velocity fields nor three-dimensional structure of the real 
space. In this way we can embed the model in a one-dimensional real space 
and simplify enormously the structure of the nonlinear terms. 
The most important advantage is that we can increase the total number of 
resolved scales (and therefore we can reach high Reynolds numbers), paying 
the price of having a model which is not exactly derivable from the original 
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NS equations. Similarities and differences between our work and the works of 

Grossmann and coworkers can be found in the conclusions. 

The main result of this paper is to present for the first time a one- dimensional 

dynamical system sharing some of the most important properties with a real 

three-dimensional turbulent flow. In particular, we analyze structure functions 

intermittency and the energy dissipation spatial distribution. RKH is found 

to be remarkably well satisfied. These models are the natural ground where 

testing new developments of large-dimensional dynamical-systems theory and 

new approaches to turbulent (i.e., multi-scale systems) flows. 

As we discuss in the following, numerical simulations need state-of-the-art 

multi-processor computers. 

The paper is organised as follows. In Sec. 2, we discuss how to jump from zero- 
dimensional shell models to one-dimensional spatially-resolved tree models; 
the tree model is introduced in Sec. 3; in Sec. 4, numerical results on structure 
functions intermittency, viscous scales dynamics and energy dissipation field 
statistics are presented; conclusions follow in Sec. 5. 



2 From chains to trees 

In order to understand turbulent energy transfer dynamics and related inter- 
mittent effects, dynamical deterministic models have been proposed. Among 
these models, shell models have recently attracted the interests of many re- 
searchers (see [[OJ for a popular introduction). The basic idea in such models 
is to retain only a few variables as representative of an entire shell of Fourier 
amplitudes. The nonlinear structure of NS equations is preserved, but all 
three-dimensional properties of the original embedding space are lost. 
The most popular shell model is the Gledzer-Ohkitani-Yamada (GOY) model 
(0-0)- 



The GOY model can be seen as a severe truncation of the NS equations. Dy- 
namical variables are described by a single complex mode u n for each shell of 
wave numbers k such as k n < k < k n+ i, with k n = A n and A being an arbi- 
trary scale parameter (A > 1), usually taken equal to 2. In the equations, only 
quasi-local couplings between nearest and next to nearest neighbor shells are 
kept. 

Recently, a new class of shell models based upon the helical decomposition 
of NS equations [I7j has been suggested [0] and studied fl^| . In this way, 



it is possible to obtain a second non-positive defined invariant closer to the 
NS helicity. The models, which have now two complex variables per shell (u+ 
and u~, transporting positive and negative helicity respectively) are simple 
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generalizations of the GOY model. From now on we will concentrate only 
on one of them (see below), which has been previously shown to share many 
properties with true turbulent three-dimensional flows [TS| . 



In particular, in [|19| it has been shown that this model has the same degree 
of intermittency found experimentally in fully developed turbulence, if the set 
of free parameters is chosen in order to conserve energy and helicity in the 
inviscid and unforced limit. 

Shell models can be thought of as field problems in zero spacial dimension: 
their obvious limitation consists in treating all degrees of freedom in a Fourier 
shell at once; the physical object they describe are coherent planar waves, filling 
the whole volume of the fluid. 

Real turbulence consists of localized eddies of all sizes that interact, merge and 
subdivide locally: the physical picture is that of a large eddy which decays into 
smaller eddies. The number of degrees of freedom in such a field problem in d 
dimensions grows with the wave number as N(k) ~ k d (d = in shell models). 
The first step in reproducing this kind of hierarchical structure is to transform 
a chain-model into a iree-model with d — 1. This is achieved by letting grow 
the number of degrees of freedom with the shell index n as 2". 
As in the original shell models, this tree model must be in some sense rem- 
iniscent of the NS equations. It can be regarded as describing the evolution 
of the coefficients of an orthonormal wavelets expansion of a one- dimensional 
projection of the velocity field v (x, t): 

v(x,t) = J2vn,j(t)i) nJ (x). (5) 

n,j 

Here ifj n j(x) are a complete orthonormal set of wavelets generated from an 
analyzing wavelet ipo$(x) by discrete translations and dilations: 

^ n ,(x) = 2 n ^ 0fi (2 n x-j). (6) 

In principle, it is possible to plug the three-dimensional generalization of (||) in 
the NS equations and deriving an exact dynamical evolution for the wavelets 
coefficients [[2C]]. For the sake of simplicity, we will be interested to some ap- 
proximated one- dimensional truncation of the wavelets resolved dynamics. In 
particular, we can construct dynamical equations of the type of shell models 
for our tree variables, viewed as the analogues of the expansion coefficients 
v n j. Similar wavelets decomposition, but with purely stochastic coefficients 
v n> j have been used for defining synthetic multi-afline signals |f2l |. 



In fig.l, we pictorially show our tree structure, covering the one-dimensional 
interval [0, At] ■ Each dynamical variable v n j is represented by a box of length 
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l n = 2~ n , occupying the region Aj-(n) ranging from (J — l)l n to jl n . At each 
scale n there are 2™" 1 boxes, covering a total length = 2 n ~ 1 l n = 1/2. 
For the sake of convenience we define the tree model in terms of density vari- 
ables, u n j (depicted as balls in figs. 2, 3), which would correspond to u n j = 
2 n / 2 v n j in a wavelets expansion. 

In this notation, |w n j| 2 represents the energy density in a flow structure of 
length l n = 2~ n and spatially labeled by the index j. 



3 The tree model 



We have chosen the helical shell model studied in [19] as the starting point for 
the construction of our hierarchical structure. The dynamical equations of this 
shell model are the following: 



W+ = ik n (a U+ +1 U n+2 + b U+_ x U n+l + C M n _ 2 Vl+)* - v K u n + S n,n F+ i ( 7 ) 

and the same holds, with all helicities reversed, for u~. Here, n = 1,...,N, 
where N is the total number of shells, v is the viscosity, F + the external 
forcing acting on a large-scale shell Uq and a,b,c are three parameters, which 
are determined by imposing conservation of energy and helicity in the inviscid 
and unforced limit: 



dE/dt = d/dt (£„ | M +| 2 + K| 2 ) = 0, 

dH/dt = d/dt (J2nkn(K\ 2 - K\ 2 )) = 0. 1 } 

The statistical properties of this model have been studied in [ l9[ : the system 
turned out to have an intermittent energy transfer very similar to what one 
can find in the original NS equations. 

Hereafter we will fix the intershell separation A = 2. For this particular (and 

standard) choice, we must fix: a = 1, b = —5/12 and c = —1/24. 

We introduce a spatial degree of freedom in the system by using the notation 

u n j to indicate the complex helical variable on scale n and spatial position 

labeled by the index j. For a given shell n, the index j can vary from 1 to 

2 n ~ 1 . 

In this tree structure, each variable u n j continues to interact with the nearest 
and next nearest levels, as in equation (0); however, a variety of possibilities 
is now opened by the presence of many horizontal degrees of freedom localized 
on each shell. 

The simplest choice is depicted in fig. 2, where a portion of the tree structure is 
showed and the evolving in time variable, u n j, is represented by a black ball. 



6 



In the figure, solid lines connect interacting balls (variables); in particular, sev- 
eral solid lines connecting one ball on a level to a group of balls, all belonging 
to another level, mean a sum over the corresponding variables in the equation. 
For simplicity, we will omit to write down explicitly these sums in the equa- 
tions, avoiding to specify the second (spatial) index j of the variables when is 
univocally determined from figures or in case of the implicit above mentioned 
sums. Thus, by writing the dynamical tree equations as follows: 

K,j = ik n {a/4: u+ +1 u~ +2 +b/2 u^u'^+c u^ 2 u~^ l )*-vk 2 n u+ j +5 riino F + , (9) 
we mean, for example, in the first term of the RHS of (13) (see fig. 2a): 

a /4 [ U n+l,2j-l( U n+2Aj-3 + U n+2,ij-2) + U n+l,2j( U n+2Aj-l + M n+2,4j)]- (10) 

Here again the same holds, with all helicities reversed, for u^j- The numerical 
values of a, b and c are the same of the original helical shell. In the unforced 
and inviscid limit our system conserves the total energy and helicity, namely: 

dE/dt = d/dt (Enj 2"1</ + K/) =0, 

dH/dt = d/dt (£ nJ 2- n k n \u+j\ 2 - K,/) = 0. { ' 

Let us remind that our tree variables can be roughly viewed as the one- 
dimensional model counterparts of the amplitudes, u n j = 2 n ^ 2 v n j, in a wavelets 
expansion of the velocity field (see eq.(H)). 

This analogy must be carried only on qualitative grounds, in particular if one 
considers that, in the wavelets basis, the Laplacian of the viscous term is not 
diagonal and that we have drastically reduced the possible interactions among 
variables. However, in a wavelets expansion, we expect the nonlinear and dis- 
sipation terms to be large only if the involved scales are of similar sizes: this 
corresponds to the physical property that eddies in a fluid mainly transfer en- 
ergy to eddies of similar scale. Also interaction among the position indices can 
be expected to be localized, with the precise degree of localization depending 
on the details of the wavelets basis. 

Besides the simple implementation (||), one can build up more complex equa- 
tions, simply turning on new interactions among variables (always conserving 
energy and helicity). The guide line in restricting the possible choices can be 
phenomenologically motivated by requiring a certain degree of locality, both in 
Fourier and real spaces. 

As an example, an enlarged set of interactions is depicted in fig.3: in particu- 
lar two contributes reminiscent of the Desnyansky-Novikov (DN) shell model 
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p2| have been added (see figs. 3d, 3e). Structurally new interactions are those 
depicted in fig.3f, where only horizontal couplings are considered. 
By putting all together we have: 



Mnj = ik n [a/A U^ +l U n+ 2 + b/2 Un-l^n+l + c u n-2 u n-l J s~ 

d {~ u n+l u n+l ~ U n+l U n+l + e l U n U n~l + e 2 U n U n-l)^~ /-j^ 
f ( U n,j+l U n,j+2 ~ U n,j-l U n,j+l ~ U n,j-l U n,j+l + U n,j-2 U n,j-l)] 

-vk 2 u + +5 F + 

lyrXj n a nj ' u n,n 1 j 

where the spatial index is shown only in the terms corresponding to the hor- 
izontal interactions of fig.3f, i.e., in the third row of eq.fli~2|). In considering 
these particular terms, the tree must be viewed as leaving on the surface of a 
cone, with the last position j on each level connected to the first one on the 
same level. 

All coefficients in the nonlinear terms of ([12]) are chosen in order to conserve 
the total energy and helicity ([□]) of the system. Each of the three interactions 
types (GOY, DN and horizontal) conserves separately E and H: this leaves 
some freedom only in the choice of the two relative weights d and / in (|T2j). 
We fixed d = f = 1; while the coefficients a,b,c are the same as before and 
ei = 3/4 and e 2 = 1/4. 

In our simulations we have considered both the simple version of the hierarchi- 
cal system described by @ and by fig. 2, which we will call version A, and the 
more complete described by (|12"D and by fig. 3, denoted hereafter as version B. 



4 Numerical results 

4.1 Numerical implementation 

In our simulations we have considered a total number of levels N = 16. The 
total number of sites forming the tree is then Nt = 2 N — 1 = 65535, each one 
described in terms of two complex variables. 

When dealing with this tree structure, the computer effort needed for numerical 
simulations increases enormously with respect to the original chain model: in 
order to collect a reliable statistics a computing time O(10 4 ) longer is now 
needed. 

For this reason, we have implemented our numerical experiments on the APE- 
100 machine [23| . 



APE-100 is a Single Instruction Multiple Data (SIMD) Parallel Architecture, 
based on a three dimensional cubic mesh of nodes with periodic boundary 
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conditions, each node being connected to its 6 neighbors. The particular version 
we used is a 512-node configuration, with a floating point performance of 50 
Mflops peak speed on a single processor. Each processing node contains a 
floating point 32 bit real arithmetics processor and a 4 Mbyte local memory. 
We used the parallelism of this machine in the simplest way: its 512 proces- 
sors integrate simultaneously the same equations, starting from different initial 
conditions. A good statistics is then obtained by time-averaging on each single 
processor and ensemble- averaging over all processors. 

In both cases A and B, we put the following parameters values into the dy- 
namical equations: 

N = 16, 

v = 2.0 10~ 5 , 

F ± = (5.0 10~ 3 ,5.0 10~ 3 ), (13) 
n = 1, 

k = 6.25 10~ 2 . 

We used the Runge-Kutta fourth-order method in integrating the system, with 
a time step varying from 10~ 3 to 5.0- 10~ 4 time units, depending on case A or B. 
In both cases we integrated the model for a total number of steps of the order 
of 10 5 on each processor of the parallel machine; considering all processors, we 
were able to collect a total statistics of 500 eddy-turnover times, using few days 
of computing time. 

4.2 Syncronization 

Before considering the statistical properties of the system, let us mention a 
dynamical effect which emerged immediately in numerical integrations, that 
is, synchronization on the last (dissipative) levels. 

Even starting to integrate from an initial configuration having different values 
on different sites of each tree level, dynamical variables on nearby locations in 
the last shells (i.e., in the dissipative range) synchronize after relatively few 
time steps. 

The mechanism underlying this phenomenon can be easily understood if one 
considers what follows: at the viscous cutoff (which corresponds in our case 
to rid ~ 11), the ratio between variables belonging to two successive scales 
l^nl/Kn-il undergoes an abrupt decrease, when passing from the inertial range 
(n < n d ) to the viscous range {n > rid). Then, in the dynamical equation for 
n = rid, interactions involving the upper levels become largely predominant 
with respect to the others: this naturally leads to a similar behaviour of vari- 
ables sitting on two close locations, whose dynamics turns out to be governed 
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by essentially the same equation. This synchronization is then transmitted and 
spread out into lower levels, in which nearby variables become equal in groups 
of 4, 8, 16, .... and so on, as the shell index n increases. 

In both cases A and B, this synchronization does not affect the inertial shells. 
In the version A of the model this process starts already from n* = 11. In 
case B the situation changes: the horizontal interactions are quite effective in 
breaking the synchronization on two close sites nearby the dissipative cutoff 
and synchronization effects are thus shifted towards the very last levels (n =15 
and 16, in our case). 

4.3 The energy dissipation field 

The first step in constructing the energy dissipation field of our model is to 
consider the following function: 

rin,j = fA£(|</ + K/), (14) 

which represents the energy dissipation density in the structure covering the 
region A 3 -(n) of length 2~ n , centered in the spatial site labeled by j. These 
structures are represented by boxes in fig. 1. 

The total energy dissipation density, e = (1/A^) J* At e(x) dx, where At is the 
total space length, is, by definition, the sum of all these contributes (sum over 
boxes at all scales in fig.l): 

e = £2-X,r (15) 

n,j 

On the other hand, in order to study the scaling properties of the energy 
dissipation field, one has to disentangle in e the contributions coming from the 
coarse-grained energy dissipation field e r , as defined in Sec.l, eq.fl2|). In our 
formulation, we can then rewrite: 




where the last expression is independent of n and the e nj 's are the coarse- 
grained energy dissipation densities, obtained as averages over spacial regions 
of length 2~ n . Note that the average density e n .j over Aj(n) does not coincide 
simply with the density r\ n j of the structure living in Aj(n), namely: 

£n,j Vn,j ^ ] Vm,k(m) ^ ] ^ Vm,k(m) ^7(m) • (1^) 
m<n m>n 



10 



Here, in the second (third) term of the RHS we take into account density 
contributions coming from larger (smaller) scale structures (as an example, all 
regions contributing to the definition of e n j are represented as shadowed boxes 
in fig.l). The index k(m) in the second term of RHS labels the location of larger 
scale structures containing the region Aj(n) under consideration (shadowed 
boxes with m < n in fig.l). In the third term, an average is performed over 
k(m) G I{m), where I(m) labels the set of structures contained in A 3 -(n), for 
any m > n (in fig.l, I(m) labels the two boxes at n+1, the four boxes at n + 2, 
and so on). 

The best spatially resolved energy dissipation field is for n = N: 

e N,j = X] Vm,k(m)'i J = l>—> 2 _1 - ( 18 ) 

m<N 

In fig. 4, the instantaneous values assumed by e^j in the N?/2 = 32768 loca- 
tions of the last level are showed. The chaotic, intermittent character of this 
spatial signal is evident. 



4.4 Scaling laws and the Refined Kolmogorov Hypothe- 
sis 

Performing long-time numerical integrations of eqs.(|9|JT2]), we have studied the 
statistical properties of both versions A and B of our tree model. We have 
investigated the scaling properties of velocity field structure functions, S p (n), 
and of coarse-grained energy dissipation moments, D p (n). 
These moments have been evaluated by time, space and ensemble-averaging 
our variables over the total integration time, the spatial locations and the 
processors of the APE machine: 



p 



S p (n) = J2 m ,t,j MT2" ( V \ U n,j{pi m )\ 2 + \ u n,j{t) m )\' 

(19) 

Dp(jl) = X)m,t,j MT2 n ^n,j{^^ m \ 

where m is the processor index (varying from 1 to M — 512) and t is the time 
step index. 

These moments are expected to follow, in the inertial range, the scaling laws 
(see eqs.(0,|)): 

S (n) = 

D p (n) =K^\ [Z[) > 
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with the two scaling exponents related together by the RKH relation (B), even- 
tually. 

In fig. 5, the log-log plot of the sixth order structure function S e (n) against 
k n is presented for the two versions A and B: the two different slopes in the 
inertial range indicate a degree of intermittency strongly dependent on how 
the variables of the tree interact. 

Let us notice that in both versions A and B, the interaction range in Fourier 
space is the same (a typical eddy at scale n interacts always with the nearest 
and nextnearest scales n — 1, n + 1, n — 2 and n + 2); nevertheless, in the two 
cases the interaction range in physical space is different and a different num- 
ber of connections between scales is considered. This seems to be of primary 
importance in determining the intermittency degree of the system. 
In order to obtain a quantitative measure of the structure functions scaling 
exponents, we performed a fit over the inertial range shells, using Extended 
Self Similarity (ESS) ||: in fig. 6, the logarithm of the sixth order moment 
5*6 (n) against the logarithm of the third order moment S 3 (n) is plotted, again 
for both versions A and B. The resulting exponents, C(p)/C(3), obtained for all 
the moments p = 1, 8, are reported in table 1 and compared with the K41 
nonintermittent prediction p/3. 

We found a similar scaling behaviour for the energy dissipation moments: as 
an example, we show in fig. 7 the log- log plot of -D 2 (n) against k n ; again, the 
reduced steepness in case B indicates a less intermittent behaviour, with r(p) 
closer to the zero value. 

The subsequent step was to test the RKH for our tree model: to improve the 
fitting procedure, we applied the ESS method to the RKH relation [fJJ: 

S p {n) = (S 3 (n)f/ 3 D p/3 (n). (21) 

Plotting the LHS of this equation against the RHS in a log-log plot, one can 
directly test the validity of relation (^) in the inertial region: in fig. 8 three 
cases, with p — 4, 5, 6, are reported, for the version A; the same, but for 
version B, are in fig. 9. 

In both versions A and B and for all the moments considered, our fits confirm 
that the RKH is well satisfied, being all the slopes in the inertial region equal 
to 1 within a few percent. 

The validity of the RKH confirms that dissipative quantities have non trivial 
influences on inertial range dynamics (and viceversa). The tree- like structure 
imposed on the velocity fluctuations does not necessary implies that the energy 
dissipation can be described in terms of fragmentation processes. In order 
to test the scale-organization of the energy structures, ultrametric-sensitive 
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observables should be studied. 



5 Conclusions 

A new class of dynamical models in one spatial dimension which show spatial 
and temporal chaos (fully developed turbulence) has been introduced and nu- 
merically studied. 

The model originates from a wavelets-like decomposition of a one-dimensional 
cut of a turbulent velocity field. Mainly local interactions in space and in scale 
are retained. 

Structure functions scaling shows typical multifractal behaviour of a real turbu- 
lent flow. Also spatial fluctuations of the energy dissipation are in qualitative 
and quantitative agreement with the corresponding measured quantities in a 
real fluid. Refined Kolmogorov Hypothesis linking the statistics of structure 
functions to the coarse-grained energy dissipation is found to hold within a few 
per cent. 

This kind of models open the possibility to investigate numerically many open 
questions of turbulence. Among them, work is in progress for detecting possi- 
ble ultrametric structure in the energy cascade mechanism |24 . 



The presence of a real space allows also to test in a quantitative way some of 
the most popular eddy-viscosity models used for simulating small-scale activity 
in Large-Eddy-Simulations. 

Our results show that intermittency strongly depends on the degree of tree 
connectiveness. For example, by passing from version A to version B one 
observes a decreasing in the deviations 5((p) = \((p) — p/3| from Kolmogorov 
scaling (from ~ 30% in case A to ~ 10% in case B): this is probably due to the 
presence of horizontal couplings in version B which allow a better efficiency in 
the energy exchanges. Increasing the number of triads-couplings should enlarge 
the number of possible downward paths followed by the energy cascade. 
This tendency toward a less intermittent regime, by increasing the number of 
triad interactions, may seem in contrast with the observed intermittency in the 
original Navier-Stokes equations (which have all possible interactions switched 
on). This contradiction is only apparent: divergenceless character of the orig- 
inal NS field, added to complicated phase-coherence effects, can very easily 
introduce different dynamical weights in the possible triad interactions leading 
to a situation where only a few of them govern the global dynamical evolution. 
For example, Grossmann and coworkers showed, by performing suitable trun- 
cation of NS equations, that intermittency depends on the typical degree of 
locality in Fourier space of the survived triad interactions; very similar results 
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have also been found in shell models at varying the inter-shell ratio A [f25] ; the 
question of finding universal observables in hydrodynamics-like intermittent 
systems remains of primary interest PS] , fZ7j ]. 
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FIGURE CAPTIONS 



Figure 1: A picture of the hierarchical system, covering the one-dimen- 
sional interval [0,At]. 

Figure 2: Pictorial representation of nonlinear interactions of eq.(||). 

Figure 3: Pictorial representation of nonlinear interactions of eq. (jl2"l) . 

Figure 4: Instantaneous configuration of the coarse-grained energy dissi- 
pation density field, ejvj, over the last level sites (version B). 

Figure 5: Log- log plot of the sixth order velocity field structure function, 
Se(n), against the wave number k n , for versions A (white circles) and B 
(black diamonds). 

Figure 6: Log- log plot of the sixth order velocity field structure function, 
Se(n), against the third order velocity field structure function, S^n), 
for versions A (white circles) and B (black diamonds). Straight lines 
correspond to a linear fit in the inertial range. The corresponding slopes 
are reported in table 1. 

Figure 7: Log- log plot of the second moment, D 2 (n), of energy dissipation 
density against the wave number k n , for versions A (white circles) and B 
(black diamonds). 

Figure 8: Log- log plot of the velocity field structure functions, S p (n) 
(p = 4, 5, 6), against ((S a (n)) p ' 3 D p / 3 (n)) in the inertial range, for version 
A. Straight lines correspond to a linear fit: all cases are compatible with 
the RKH. 

Figure 9: The same as in fig. 8, but for version B. 

TABLE CAPTION 

Table 1: The scaling exponents ratios, C(f0/C(3) (p — 1, 8), resulting 
from the ESS fitting procedure. The values obtained for versions A and 
B are reported in the second and third column, respectively. In the last 
column there are the K41 values. 
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V 


A 


B 


KAl 


1 


0.41 ±0.01 


0.348 ± 0.005 


0.333 


2 


0.74 ±0.01 


0.682 ±0.005 


0.667 


3 


1 


1 


1 


4 


1.21 ±0.01 


1.303 ±0.006 


1.333 


5 


1.36 ±0.02 


1.59 ±0.01 


1.667 


6 


1.48 ±0.03 


1.86 ±0.02 


2 


7 


1.55 ±0.05 


2.12 ±0.03 


2.333 


8 


1.60 ±0.07 


2.35 ±0.03 


2.667 
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